Global well-posedness of 2D compressible 
Navier-Stokes equations with large data and vacuum 



Quansen Jiu, 1,3 * Yi Wang 2,3 ^ and Zhouping Xin 3 ^ 



1 School of Mathematical Sciences, Capital Normal University, Beijing 100048, P. R. China 

2 Institute of Applied Mathematics, AMSS, and Hua Loo-Keng Key Laboratory of Mathe- 
matics, CAS, Beijing 100190, P. R. China 

3 The Institute of Mathematical Sciences, Chinese University of HongKong, HongKong 

Abstract: In this paper, we study the global well-posedness of the 2D compressible Navier- 
Stokes equations with large initial data and vacuum. It is proved that if the shear viscosity 
/i is a positive constant and the bulk viscosity A is the power function of the density, that is, 
A(p) = with /3 > 3, then the 2D compressible Navier-Stokes equations with the periodic 
boundary conditions on the torus T 2 admit a unique global classical solution (p, u) which may 
contain vacuums in an open set of T 2 . Note that the initial data can be arbitrarily large to 
contain vacuum states. 
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1 Introduction 

In this paper, we consider the following compressible and isentropic Navier-Stokes equations 
with density-dependent viscosities 

dtp + div(pu) = 0, 

d t {pu) + div(pu <g) u) + VP(p) = pAu + V((p + A(p))divu), x E T 2 ,t > 0, 

where p(t,x) > 0, u(t,x) = (ui,U2)(t,x) represent the density and the velocity of the fluid, 
respectively. And T 2 is the 2-dimensional torus [0, 1] x [0, 1] and t £ [0, T] for any fixed T > 0. 
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We denote the right hand side of ([l.lP o by 



C p u = pAu + V((/i + A(p))divu). 



Here, it is assumed that 



H = const. > 0, 



Hp) = p< 



13 > 3 



(1.2) 



such that the operator C p is strictly elliptic. 
Let the pressure function be given by 



P(p) 



(1.3) 



where 7 > 1 denotes the adiabatic exponent and A > is the constant. Without loss of 
generality, A is normalized to be 1. The initial values are given by 



Here the periodic boundary conditions on the unit torus T 2 on (p,u)(t,x) are imposed to the 
system (jl.ip . This model problem, (|l.ip - (|1.4p . was first proposed by Vaigant-Kazhikhov in 
|51j where they showed the well-posedness of the classical solution to this problem provided the 
initial density is uniformly away from vacuum. In this paper, we study the global well-posedness 
of the classical solution to this problem (|l.ip - (|1.4p with general nonnegative initial densities. 

There are extensive studies on global well-posedness of the compressible Navier-Stokes equa- 
tions in the case that both the shear and the bulk viscosity are positive constants satisfying 
the physical restrictions. In particular, the one-dimensional theory is rather satisfactory, see 
[20\ [38] [33l [M] and the references therein. In multi-dimensional case, the local well-posedness 
theory of classical solutions to both initial-value and initial-boundary-value problems was es- 
tablished by Nash [H], Itaya [26] and Tani [50] in the absence of vacuum. The short time 
well-posedness of either strong or classical solutions containing vacuum was studied recently by 
Cho-Kim [8] and Luo[50] in 3D and 2D case, respectively. In particular, Cho-Kim [8] obtained 
the short existence and uniqueness of the classical solution to the Cauchy problem for the isen- 
tropic CNS with general nonnegative initial density under the assumption that the initial data 
satisfies a natural compatibility condition [8] . One of the fundamental questions is whether these 
local (in time) solutions can be extended globally in time. The first pioneering work along this 
line is the well-known theory of Matsumura-Nishida [41] . where they obtained a unique global 
classical solution to the CNS in .fP(M 3 ) (s > 3) for initial data close to its far field state which 
is a non-vacuum equilibrium state, and furthermore, the solution behaves diffusively toward the 
far field state. The proof in [41] consists of elaborate energy estimates based on the dissipa- 
tive structure of the CNS and spectrum analysis for the linearized of CNS at the non-vacuum 
far field state. This theory has been generalized to data with discontinuities by Hoff [18] and 
data in Besov spaces by Danchin in [9]. It should be noted that this theory [41[ [TB] [9] requires 
that the solution has small oscillations from the uniform non-vacuum far field state so that 
the density is strictly away from the vacuum uniformly in time. A natural and important long 
standing open problem is whether a similar theory holds for the initial data containing vacu- 
ums. In this direction, the major breakthrough is due to P. L. Lions [37] . where he obtained 
the existence of a renormalized weak solution with finite energy and large initial data which 
can contain vacuums for the isentropic CNS when the exponent 7 is suitably large, see also the 
refinements and generalizations in [15\ 129]. However, little is known on the structure, regularity, 
and uniqueness of such a weak solution except the partial regularity estimates for 2-dimensional 



(p,u)(t = 0,x) = (p ,u ){x). 



(1.4) 
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periodic problems in Desjardins [10] where a stronger estimate is obtained under the assump- 
tion of uniform boundedness of the density. Recently, under some additional assumptions on 
the viscosity coefficients, and the far fields state is a non-vacuum state, Hoff |18[ [TU] obtained 
a new type of global weak solution with small total energy for the isentropic CNS, which have 
extra structure and regularity information (such as Lagrangian structure in the non-vacuum 
region) compared with the renormalized weak solutions in [37 \ I15 [ [29]. However, the uniqueness 
and regularity of those weak solutions whose existence has been proved in |37[ \TE\ [29] remain 
completely open in general. By the weak-strong uniqueness of P. L. Lions [37], this is equivalent 
to the problem of global (in time) well-posedness of classical solution in the presence of vacuum. 
It should be pointed out that this important question is a very difficult and subtle issue since, in 
general, one would not expect a positive answer to this question due to the finite time blow-up 
results of Xin in [52], where it is shown that in the case that the initial density has compact 
support, any smooth solution to the Cauchy problem of the CNS without heat conduction blows 
up in finite time for any space dimension, see also the recent generalizations to the case for 
non-compact but rapidly decreasing (at far fields) initial density [36]. The mechanism for such 
a blow-up has also been investigated recently and various blow-up criterion have been derived 
in [El [IH 1221 [23 123 SSI SS] - More recently, Huang-Li-Xin[24J proved the global well-posedness 
of classical solutions with small energy but large oscillations which can contain vacuums to 3D 
isentropic compressible Navier-Stokes equations. See also the recent generalization to 3D full 
compressible Navier-Stokes equations [21] . the isentropic Navier-Stokes equations with potential 
forces [35], and ID or spherically symmetric isentropic Navier-Stokes equations with large initial 
data HH Q2J. 

The case that the viscosity coefficients depend on the density and vanish at the vacuum has 
received a lot attention recently, see [21 13 SI 13 [SI IIZ1 12Z1 123 123 ffl 1^ 1321 SSI S21 1^1 SI1 
[53| l54l [55] and the references therein. Liu, Xin and Yang first proposed in [39] some models of 
the compressible Navier-Stokes equations with density-dependent viscosities to investigate the 
dynamics of the vacuum. On the other hand, when deriving by Chapman-Enskog expansions 
from the Boltzmann equation, the viscosity of the compressible Navier-Stokes equations depends 
on the temperature and thus on the density for isentropic flows. Also, the viscous Saint- Venant 
system for the shallow water, derived from the incompressible Navier-Stokes equations with a 
moving free surface, is expressed exactly as in (11. ip iV = 2, fi = p, A = 0, and 7 = 2 (see 
[16]). For the special case, (|1.2p . the global well-posedness result of Vaigant-Kazhikhov [51] 
is the first important surprising result for general large initial data with the only constraint 
that it is initially away from vacuum. However, in the presence of vacuum, there appear new 
mathematical challenges in dealing with such systems. In particular, these systems become 
highly degenerate. The velocity cannot even be defined in the presence of vacuum and hence 
it is difficult to get uniform estimates for the velocity near vacuum. Substantial achievements 
have been made for the one-dimensional case, such as both short time and long time existence 
and uniqueness for the problem of a compact of viscous fluid expands into vacuum with either 
stress free condition or continuity condition have been established with p = p a for suitable a, 
see [23 E3 (S3 El] etc. Li-Li-Xin [36] recently proved the global existence of weak solutions 
to the initial-boundary value problem for such a system on a finite internal with general initial 
data which may contain vacuum and discovered the phenomena that all the vacuum states must 
vanish in finite time and any smooth solution blows up near the time of vacuum vanishing which 
are in sharp contrast to the case of constant viscosity coefficients, which have been extended 
to the Cauchy problem on M 1 for arbitrary initial data with a uniform non-vacuum state at 
far fields by Jiu-Xin [32]. In the case that a basic nonlinear wave pattern is the rarefaction 
wave, whose nonlinear asymptotic stability has been proved in [301 [3"T] for the one-dimensional 



4 



isentropic CNS system with density-dependent viscosity in the framework of weak solutions 
even the rarefaction wave is connecting to the vacuum[38]. Note also that in the case that the 
initial data is strictly away from vacuum, Mellet and Vasseur has obtained the existence and 
uniqueness of global strong solution to the one-dimensional Cauchy problem [33]. However, the 
progress is very limited for multi-dimensional problems. Even the short time well-posedness 
of classical solutions has not been established for such a system in the presence of vacuum. 
The global existence of general weak solutions to the compressible Navier-Stokes equations with 
density-dependent viscosities or the viscous Saint- Venant system for the shallow water model 
in the multi-dimensional case remains open, and one can refer to [4], [5J, [17], [42] for recent 
developments along this line. Note also that Zhang-Fang [55] proved the existence of global 
weak solution with small energy to 2D Vaigant-Kazhikhov model [51] in the framework of [19] 
and presented the vanishing vacuum behavior. However, the uniqueness of this weak solution is 
open. 

In this paper, we investigate the global existence of the classical solution to 2-dimensional 
Vaigant-Kazhikhov model [51] . that is, CNS system (|l.ip - (|1.4p with periodic boundary condition 
and general nonnegative initial density. It should be noted that for the 2-dimensional problem, 
the basic reformulation of Vaigant-Kazhikhov [51] and the formulation in terms of the mate- 
rial derivative used in [18} [24] are equivalent. Following some of the key ideas developed by 
Vaigant-Kazhikhov [51], we are able to derive the uniform upper bound of the density under 
the assumptions that the initial density is nonnegative. Then we can derive the higher order 
estimates to the solution to guarantee the existence of the global classical solution. 

The main results of the present paper can be stated in the following. 
Theorem 1.1 If the initial values (po,uq)(x) satisfy that 

< ( Po (x),P(p )(x)) G W 2 > q (T 2 ) x W 2 '«(T 2 ), u (x) G H 2 (T 2 ), [ Po (x)dx > (1.5) 
for some q > 2 and the compatibility condition 

£ P0 u - VP(po) = Vp g(x) (1.6) 

with some g G L 2 (T 2 ), then there exists a unique global classical solution (p,u)(t,x) to the 
compressible Navier-Stokes equations (|l.ip - (|1.4p with 

0<p(t,x)<C, y(t, x) G [0, T] x T 2 , ( ft P(p))(^)eC([0,T];^«(T 2 )), 
u G C([0,T];H 2 (T 2 ))r\L 2 (0,T;H 3 (T 2 )), y/tu G L°°(0, T; H 3 (T 2 )), 

tu G L°°(0, T; W 3 ' g (T 2 )), u t G L 2 {0, T; H 1 ^ 2 )) (1.7) 

Vtu t G L 2 (0,T;F 2 (T 2 ))nL oo (0,T;i7 1 (T 2 )), tu t G L°°(0, T; H 2 (T 2 )), 

Viy/putt G L 2 (0,T;L 2 (T 2 )), t^pu tt G L°°(0, T; L 2 (T 2 )), tVu tt G L 2 (0, T; L 2 (T 2 )). 

Remark 1.1 From the regularity of the solution (p,u)(t,x), it can be shown that (p,u) is a 
classical solution of the system (11. ip in [0,T] x T 2 (see the details in Section 5). 

Remark 1.2 If the initial data contains vacuum, then it is natural to impose the compatibility 
(jl.6p as the case of constant viscosity coefficients in f^j. 
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Remark 1.3 In Theorem ! 1 . 11 it is not clear whether or not uu E L 2 (0, T; L 2 (T 2 )) even though 
one has the regularity tVu u E L 2 (0, T; L 2 (T 2 )). 

Remark 1.4 It is open to get the similar theory to the Cauchy problem or the Dirichlet problem 
to the 2D compressible Navier-Stokes equations (jl.ip . 

If the initial values are much more regular, based on Theorem 11.11 we can prove 
Theorem 1.2 If the initial values (pq,uq)(x) satisfy that 



and the compatibility condition (|1.6|) . then there exists a unique global classical solution (p, u)(t, x) 
to the compressible Navier-Stokes equations (|l.ip - (|1.4p satisfying all the properties listed in (|1.7p 
in Theorem with any 2 < q < oo . Furthermore, it holds that 



Remark 1.5 In fact, the conditions on the initial velocity uq can be weakened to uq E H 2 (T 2 ) 
and VA)V 3 u E L 2 (T 2 ) to get Q|> . 

Remark 1.6 In Theorem \1.2l it is not clear whether or not u E C([0, T]; H 3 (T 2 )) even though 
one has pu E C([0, T\; H 3 (T 2 )). 

Remark 1.7 It is noted that in Theorem \1.2\ the compatibility condition (jl.6p is exactly same 
as in Theorem \l.l\ 

Notations. Throughout this paper, positive generic constants are denoted by c and C, which 
are independent of 5, m and t E [0, T], without confusion, and C(-) stands for some generic 
constant(s) depending only on the quantity listed in the parenthesis. For function spaces, 
L P (T 2 ), 1 < p < oo, denote the usual Lebesgue spaces on T 2 and || • || p denotes its LP norm. 
W k 'P(T 2 ) denotes the k th order Sobolev space and H k (T 2 ) := W k > 2 {T 2 ). 

2 Preliminaries 

As in |51j . we introduce the following variables. First denote the effective viscous flux by 



< (po(x), P(po)(x)) E H 3 (T 2 ) x H 3 (T 2 ), u Q (x) E H 3 (T 2 ) 




(1.8) 



u E L 2 (0,T;H 4 (T 2 )), (p,P(p)) E C([0, T]; H 3 (T 2 )), 
pueC([0,T};H 3 (T 2 )), ^5V 3 u E C([0, T]; L 2 (T 2 )). 



(1.9) 



F = (2/i + A(p))div«-P(p), 



and the vorticity by 



w = 5 Xl n 2 - 9 X2 tti. 



Also, we define that 



H = 



-(pui Xl + F X2 ) 



L = 



-(-pu X2 + F X1 ). 
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Then the momentum equation (jl.ip n can be rewritten as 



u u + u- Vui = ^(-pu X2 +F X1 ) = L, 
u 2t + u- Vu 2 = j((JW Xl + F X2 ) = H. 

Then the effective viscous fiux F and the vorticity u solve the following system: 
cat + u ■ Vw + udivu = H Xl — L X2 , 

(§5$)* + u ■ V (fS$) + («i-i) 2 + 2u 1X2 U2 Xl + (u 2x2 ) 2 = H X2 + L X1 . 
Due to the continuity equation (j 1 . 1 j) x . it holds that 

ut + u ■ Vw + wdivu = H Xl — L X2 , 

F + u-VF- p(2p + A (p))[ J P(^ j y + (jjg^yjdivti 

+(2p + X(p))[{u lxi ) 2 + 2u lx2 u 2xi + {u 2x2 ) 2 } = (2p + X(p))(H X2 +L Xl ). 

Furthermore, the system for (H, L) can be derived as 



(2.1) 



(2.2) 



(2.3) 



pH t + pu • VH — pHdivu + u X2 ■ VF + pu Xl • Vco + p(udivu) Xl 
-{p(2p + X(p))[F( * )' + (^y]divn} a 



■•■■> 



(2.4) 



+{(2(i + X(p))[(u lxi ) 2 + 2u 1X2 u 2xi + (u 2x2 ) 2 ]} X2 
= [(2A* + X{p)){H X2 + La-Jk;, + ~ A^zi, 

pLt + pu • VL — pLdivu + • VF — jUit^ • Vo; — p(u;divu) X2 
-{ P (2p + A(p))[n^)' + (^Jldivu}^ 
+ {(2/i + X(p))[(u lxi ) 2 + 2u 1X2 u 2xi + (u^a) 2 ]}^ 
= [(2/i + KP)){H X2 + L X1 )] X1 - p(H Xl - L X2 ) X2 . 

In the following, we will utilize the above systems in different steps. Note that these systems 
are equivalent to each other for the smooth solution to the original system (jl.ip . 

Several elementary Lemmas are needed later. The first one is the Gagliardo-Nirenberg 
inequality which can be found in 



Lemma 2.1 V/i G W 1,m (T 2 ) or h G W l ' m (T 2 ) with / hdx = 0, it holds that 

Jt 2 



< clival 



a m \\h\\l~ a , 



(2.5) 



where a = (- — -)(- — — + 5) > an d ifm<2, then q is between r and 53^; ^hat is, q G [r, 53^] 
i f r < 1 G t^> r ] ^ ifm = 2, then q G [r,+oo), if m > 2, then q G [r,+oo]. 

Consequently, V/t G W 1,m (T 2 ), one /ias 



< + || V^ 



"Hi, ill— a - * 



(2.6) 
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where C is a constant which may depend on q. 
The following Lemma is the Poicare inequality. 

Lemma 2.2 V/i £ W 1,m (T 2 ) or/iG W 1 > m (T 2 ) / /idx = 0, i/ 1 < m < 2, f/ierc 

Jt 2 



2m < C(2 - m)"2 ||V/i|| m , (2.7) 

2— m 

where the positive constant C is independent of m. 

The following Lemma follows from Lemma 12.21 of which proof can be found in [51] . 
Lemma 2.3 V/i £ ' m + 1 ' (T ) urei/t m > 2 and < 77 < 1, we have 



2m < CiWhh + in^WhW^^WVhW 1 ^), (21 

m+77 



where e £ [0, ^],s = m-^i-ef i^e positive constant C is independent of m. 

3 Approximate solutions 

In this section, we construct a sequence of approximate solutions by making use of the theory 
of Vaigant-Kazhikhov [51] and derive some uniform a-priori estimates which are necessary to 
prove Theorem II .li To this end, we need a careful approximation of the initial data. 

Step 1 . Approximation of initial data: To apply the theory of Vaigant-Kazhikhov [51] , we 
approximate of the initial data in (jl.8p as follows. First, the initial density and pressure can be 
approximated as 

P 5 = P0 + S, P$ = P( Po ) + 5, (3.1) 

for any small positive constant 5 > 0. To approximate the initial velocity, we define Uq to be 
the unique solution to the following elliptic problem 

£ p ^o = Vitf + VP^9 (3.2) 

with the periodic boundary conditions on T 2 and / u^dx = j u^dx := uq. It should be noted 

that Uq is uniquely determined due to the compatibility condition (jl.6p . 
It follows from (|3.2|) that 

£po4 = -V [(AU) - A(po))div«g] + VP 5 + Vw?. (3.3) 
By the elliptic regularity, it holds that 

\\ U - ^o\\h 2 (T 2 ) 



< c 

< c 

< c 



\\M) - A(po)||oo||V(divt*g)[| 2 + ||V(A(pg) - A(po))||oo||div«g|| 2 + ||Vitf || 2 + \\^gh 

^ll n oll-ff 2 (T 2 ) + II-PoIIh^t 2 ) + IIv^IIl 00 ^ 2 )!!^!^ 

(T 2 ) 



^II w oIIh 2 (t 2 ) + 1 



(3.4) 
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where the generic positive constant C is independent of 5 > 0. 
Therefore, if 5 <C 1, then (|3.4p yields that 

||«oll/P(T2) < C (3-5) 
where the positive constant C is independent of < 5 <C 1. 

Due to the compatibility condition f 1 1 . 6 1) and (|3.2p . it holds that 

C po (4 - n ) = -V [(X(p s ) - A(p ))divng] := O 5 . (3.6) 
Therefore, by the elliptic regularity, (|3.ip and (|3.5p . one can get that 

\\ u ~ u o||//2(T2) < C||© 5 ||2 



< c 

< c 



WKpI) - A(po)||l-(t2)||V 2 ^|| 2 + ||V(A(pg) - A(po))||L»(^)||div^|| 2 
WKpq) ~ A(po)IU»(r«) + l|V(A(pg) - A(po))IU»(ra) 



(3.7) 



< C5 -> 0, as 5 -> 0. 



For the initial data (/Oq,Pq,Uq) constructed above for each fixed 5 > 0, it is proved in [51] that 
the compressible Navier-Stokes equations (jl.ip with /3 > 3 has a unique global strong solution 
(p s ,u s ) such that c§ < p s < C$ for some positive constants cs,C$ depending on S. In the 
following, we will derive the uniform bound to (p s ,u s ) with respect to 5 and then pass the limit 
5 to get the classical solution which may contain vacuum states in an open set of T 2 . It 
should be noted that in comparison with estimates presented in |51j . we will obtain uniform 
estimates with respective to the lower bound of the density such that vacuum is permitted in 
these estimates. To this end, the compatibility condition (ll.6p will be crucial. 

For simplicity of notations, we will omit the superscript 5 of (//, u s ) in the following in the 
case of no confusions. 

Step 2. Elementary energy estimates: 
Lemma 3.1 There exists a positive constant C depending on (po,uo), such that 

sup (HVHI 2 + \\pm) + f T (HVuH 2 , + |M| 2 + ||(2u + A(p))5divu|| 2 )dt < C. (3.8) 
te[o,T] J o 

Proof: Multiplying the equation (|2.ip , by pui, (i = 1,2), summing the resulting equations and 
then integrating over T 2 and using the continuity equation (|l.ip i . it holds that 

p\u\ 2 dx + J (puj 2 + (2p + \(p))(di\uf)dx + J u-VPdx = 0. 

Multiplying the continuity equation (|l.ip i by ^rj/J 7_1 and then integrating over T 2 yields that 

d r p* , f 

dx + Pdivudx = 0. 



dt J 7-1 

Therefore, combining the above two estimates and then integrating over [0, t] with respect to t, 
we obtain 

/ (\p\ u \ 2 + T~[P^ dx + J + ^ + Hp))(&vu) 2 )dxdt 



< c 



5 U H 3 {J 2 )\\ U o\\h2(J2) + ||/5ol|^-3( T 2) 



(3.9) 

< C. 
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Denote 



Then 



cj){t) = / (pu 2 + (2/x + X(p))(divu) 2 )dx, t £ [0,T]. 



(3.10) 



||V«|||(i) < C[\\cj\\ 2 2 (t) + ||div«||l(t)] < C<j>{t) e L\0,T). 

Thus the proof of Lemma 13. II is completed. □ 

Step 3. Density estimates: Applying the operator div to the momentum equation (jl.ip o. we 
have 

[div(/0u)]t + div[div(/9ii ® it)] = AF. 
Consider the following two elliptic problems: 



A£ = div(pu), J £dx = 0, 



At] = div[div(pit ® it)], 



7]dx = 0, 



both with the periodic boundary condition on the torus T 2 . 
By the elliptic estimates and Holder inequality, it holds that 

Lemma 3.2 (1) ||V£||2m < Cm||p|| 2mk ||it||2mfc) f or an V k > l,m > 1; 



(2) ||V£|| 2 _ r < C\\Jpu\\ 2 \\p\\U, for any < r < 1; 

r 

(3) \\r)hm < Cm\\p\\2mh\\u\\l mk , for any k > l,m> 1; 

A; — 1 

where C are positive constants independent ofm,k and r. 



(3.11) 



(3.12) 



(3.13) 



Proof: (1) By the elliptic estimates to the equation (|3. 12[) and then using the Holder inequality, 
we have for any k > 1, m > 1, 



||V£|| 2 m < Cm\\pu\\ 2m < Cm\\p\\2mk \\u\ 

fe-i 

Similarly, the statements (2) and (3) can be proved. 

Based on Lemmas I2.1H2.3I and Lemma 13.21 it holds that 



Irak ■ 



□ 



Lemma 3.3 (1) ||£||2m — Cm? ||V£|| 2m < Cmz \\p\\m, for any m > 2; 

m + 1 



(2) \\u\\ 2m < C 



ma llVulU + 1 



for any m > 2; 



(3) ||V£||2m < C 1712 \\p\\ 2mfc (j)(t) 2 +m||/0||2mfe 



for any k > 1, m > 1; 



(4) 1 1 ?7 1 1 2m < C m 2 k\\p\\ 2mfe (j)(t) + 2mt , for any k > l,m > 1; 

where C are positive constants independent ofm,k. 
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Proof: (1) By Lemma [2T21 and Lemma [3721 (2). it holds that 



i i I 

p\\m 

m+l 



< Cm*\\p\\?k, 

where in the last inequality one has used the elementary energy estimates (|3.9p . 

(2). From the conservative form of the compressible Navier-Stokes equations and the 
periodic boundary conditions, we have 

— J p(t,x)dx = — J pu(t,x)dx = 0, 

that is, 

J p(t,x)dx = J po(x)dx, J pu(t,x)dx = J poUo(x)dx, Vi £ [0,T]. 
By Lemma 12,2} it follows that 

i 

M 2m < \\u - u 2m + Hbm - Cm 1 * \\Vu\\ am + it , (3-14) 

m+l 

where m > 2 and -u = u(i) = ^ u(t,x)dx. 
On the other hand, we have 

| / p(u — u)cfcc| < ||/9|L||it — rZll _2_ < C||Vu||2, (3.15) 
J t- 1 



where in the last inequality we have used the elementary energy estimates (|3.9p and the Poincare 
inequality. 

Note that 



p(u — u)dx\ = | J pQUQdx — u J po(x)dx\ > \u\ J podx — \ J poUodx\. (3.16) 

Combining f)3. 15|) with (|3,16p implies that 

\Jp u dx\ + C^\ 2 _ 
J podx J p dx 

Substituting (|3.17p into (|3.14j) completes the proof of Lemma [3/3] (2). 

The assertions (3) and (4) in Lemma 13.31 are direct consequences of Lemma 13.31 (2) and 
Lemma 13.21 (1). (3), respectively. Thus the proof of Lemma 13.31 is completed. □ 

Substituting (ETT2]) and (l37T3|) into (|3TTT]) yields that 

A(f t + »7-F + J ' F(t,x)dxj =0. (3.18) 

Thus, it holds that 

& + V-F + J F(t,x)dx = 0. (3.19) 
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It follows from the definition of the effective viscous flux F that 

& - (2/x + A(p))div« + P(p) + V + J F(t, x)dx = 0. 
Then the continuity equation (jl.ip i yields that 



Define 



P V 



9{ P ) = P^±^ds = 2»\n P + 



Then we obtain the following transport equation 

(£ + 0(p))t + « • V(C + 0(p)) + P(p) + v - u ■ V£ + y F(t, x)dx = 0. 



(3.20) 

(3.21) 
(3.22) 

(3.23) 



Lemma 3.4 For any k > 1, it holds that 



sup ||p(V)IU <CkT- 
te[o,T] 



(3.24) 



Proof: Multiplying the equation ([3T23|> by />[(£ + ^p))^ 2 " 1 " 1 with m > 4 being integer, here 
and in what follows, the notation (•••)+ denotes the positive part of (• • • ), one can get that 



~ j P[{i + 0{p)) + ] 2m dx + J pP(p)[(t + 0(p))+] 



dx 



pri[(Z + 0(p)). 



l2m-l 



dx - / F{t,x)dx / p[(£ + 9{p)) 



i2m-l 



dx. 



Denote 



/(*) = { J p[(Z + e(p)) + ] 2m dx}^, t€[o, 

Now we estimate the terms on the right hand side of (|3.25p . First, 

p V i(t + 0(p)) + ] 2m - ^x] < [ pMv\[p(Z + e( P ))l m ] 2r 



dx 



< 



\p\\2^+l\\v\\2m+dP^ + 9 (p)) 



2m 1 1 2m 

+ 111 



< C\\p 



1 

I 2m 

\2m0+l 



1 



^0(t) + (m+ — )\\p\\^ 



f(t) 



2m- 1 



Next, for 



)]) and in the 1 
1 with p,q > 1, one has 



where <p(t) is defined as in (|3.10p and in the last inequality we have taken k = 



i +1 + 1 

2m/3+l ' p ' q 



I 



pu ■ V£[(£ + 9{p)) + \ 2 ™- 1 dx\ < J p^ |u||V£| [p(£ + 0(p)) 
< \\p\\^ + i\\ u h mp \N^h mq \\p^ + 0(p)) 



2m-l 

*rn\ 2m 



2m — 1 
\2mil 2m 



1 



-L ^ 3 1 1 

< c\\ 

Pll2m/3 + l ( m P) 2 II II 2 + 1 C™/) 2 k 2 \\p\\ 2mqk g>(t) 2 + m||p|| 2m 9 fc f (t) 



\2m-l 



<C|Ml^i/(*)^>V(t) + ro], 



dx 

(3.25) 
(3.26) 



(3.27) 



(3.28) 
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where in the third inequality one has chosen p = q = 2 "^g l ~ 1 and k - 



/3-r 



Then it follows that 

F(t,x)dx J p[(t + e(p)) + ] 2m ~ 1 dx\ 

< J \{2 l L + \{p))&\vu-p(p)\dx [ P ^[p{i + e{p)) 

< 

< c 

< c 



2m-l 
2ml 2m 



dx 



{J {2p + \{p)){d\wu) 2 dx) ± 2{J {2p + \(p))dx)i + / P{p)dx ||p|iri|p(£ + #(p)) 
'^t) 1 * +mhj P^dx) 1 * +l]f(t) 2m - 1 

'm l +mhp\\u + i+ i }fi t ? m ' x - 

Substituting (ET271) . (ET28]) and (pT29l) into (ET25|) yields that 

^^(/ 2m W)+ / pPipM + Oip^+r^dx 
2m dt J 



2m I 
+ I 



< 



^iipii^+i/^) 2 ™" 1 +m\+c U(t)a + 0(*)5|| P ii^9 +1 + ij /(*) 



\2m-l 



Then it holds that 

2m/3+l + { m ' 2( l ) (t) + ||p||2 



i 

2m 

2m/3+l 



Integrating the above inequality over [0,t] gives that 



1+^ HT^\\p\\Lp +1 (T)dT + ^ 



(mVfr) + ™)l|p|l2m#fi( T ) d ' 



Now we calculate the quantity 

/(o) = (|p^[(^ + ^)) + ] 2m ^)^- 

By Lemma 13.21 (1) with t = 0, we can easily get 

||£olk°° ^ C- 

Furthermore, by the definition of 0(p s ) = 2/ilnpg + ^{{plf - 1), we have 

^ + 0(pg)^-oo, as 4^0 + . 
Thus there exists a positive constant a, such that if < p^ < a, then 

(£o + %o))+ = 0. 

Now one has 



no) = / + / k^ + %*)) 2 ™^ 

L V J [0<po<cr] J [<r<Po<M] ' 

p 5 M + 0{pl)f+dx]^ <C(a,M), 

r<P 5 <M] J 

where the positive constant C(a,M) is independent of 5 and m. 



i 

2/n 



(3.29) 



(3.30) 



(3.31) 



(3.32) 



(3.33) 
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It follows from (137321) and (13331 that 



f(t) < C 



1 + 



o 



(r)Hp\\L^ 



1+7 



(r)dr + / (m>(r) + m)||p|| 2 ^'? 1 (r)dr 



(3.34) 



Set fii(t) = {x £ T 2 |p(t,a;) > 2} and 2 (i) = {x £ ^(rj)^, x) + 0(p)(t, x) > 0}. Then one has 
llp||£n*-i(*) = 



< 



2m/3+l 



fil(t) 



Hi(t) 

c( [ P \e(p) + i-i\ 2m dx+ f 

y JVL 2 (t) JO 



fil(t) 



2mj„ 2,n/3+l 



T 2 \fii(t) 

/3 



Oi(t)\n 2 (t) 



2m I \ 2m^+l 



(3.35) 



< c 



n 2 (t) 

\2m 



p{9(p) + ty m dx+ / P \z\ 2m dx + 



<c(f(ty m + / p|<ei 2m (ix 



2mjJ 2m/3 + l 



i9 



n 2 (t) 



2m > \ 2m,8 + l 



ni(t)\n 2 (t) 



+ c < c 



~f{t)+(j^p\t\ 2m dxy 



+i 



Note that 



2m I \ 2m/3 + l 



T 2 



s 







< \\p\\^A\tr\\zs& = \\p\\zs+i\\t\\ 

JL_ r 1 



< 



I 2m/3 + l 



1m.fi 

2m@+l II All 2m/3 + l 

2m/3+lll?H 2m +l 

2m/3 



/3(m + l) 



2 ^ +1 <Cml lipids, 



(3.36) 



Then one can get 



£(m+l) 

iJ^+i(*)^^i + /(*) + ^iipIISS+i(*) 



< 



1, 



I2771/M 



(t) + C l + /(t) + m '"(2/3-1) 



(3.37) 



Thus it holds that 



< C 



m 2/3-1 _|_ 



< C 



m 2 /3-l + 



W+i( r ) dr+ / (^(r) + rn)\\p\\ 2 l^ 1 (r)dr 



Applying Gronwall's inequality yields that 

\\p\\Lp+i(f)<c 



m 2/3-l 



.1+3=; 



Denote 

Then it holds that 



/(*) < c 

< c 



m>(T)+m)||p|| 2m 2™ 1 (r)d 
y(t) =m~3=T||p|| 3m ^ + i(*). 

/3(l-3/3) 1 /■* 1 , 1 l -i ft 1,1- 

i) _|_ m m(M) / (/>(r)y(r) i+ 2^(iT + m" 1 ^" 1 ) / y( r ) 



(3.38) 
(3.39) 



m (2/3-l)(/3 



1 + / (0(r) + l)/(r)dr 
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So applying the Gronwall's inequality again yields that 

y(t)<C, Vte[0,T], 

that is, 

\\p\\2mP+i{t)<CmT&, Vt€[0,T\. 
Equivalently, (|3.24|) holds. Thus Lemma 13.41 is proved. 
Step 4: First-order derivative estimates of the velocity. 

Lemma 3.5 There exists a positive constant C , such that 



□ 



sup / (poo 2 + 
te[o,T] J 



F 2 



2fi + X(p) 



)dx + J J p{H 2 + L 2 )dxdt < C. (3.40) 

, respectively, and 



Proof: Multiplying the equation (|2.3p i by p,oj, the equation (|2.3|) 2 by 2 t i+\( p ) 
then summing the resulted equations together, one has 



1 d 
2dt 



[pu 2 + 



■3' 



2/i + X(p) 



divu 



)dx + — [ uj 2 d\vudx / pF 2 ( — — Ydivudx 

2 J 2 J 2p + X(p) 

P(P) v. 



2p + X(p) 
J p(H 2 + L 2 )dx. 



dx —J pF(divu)( 



2p + X(p)' 



)'dx+ I F[(ui Xl ) + 2ui X2 u 2xi + (u 2x2 ) ]dx 



(3.41) 



Notice that 

\2 i o /„. \2 



(ui Xl ) + 2ui X2 u 2xi + [u 2x2 ) = (ui Xl + u 2x2 ) + 2(ui X2 u 2xi - ui Xl u 2x 



(divu) 2 + 2(u lx2 u 2xi - u lxi u 2x2 ) 



(divu) 



F + P(P) 



2p + X(p) 2p + X(p) 



+ 2(ui X2 u 2xi - ui Xl u 2 



then one has 
\jtl^ 2 + 



F 2 



2p + X(p) 



)dx 



+ J p(H 2 + L 2 )dx = ~ j uj 2 di\udx 



+ ^ j F 2 (divu) p( 



2p + X{p)' 2p + X(p) 
J 2F(u 1x2 u 2xi ~ ui Xl u 2x2 )dx. 



dx + J F(dfvu) p( 



P(P) w P(P) 



dx 



2p + X{p)' 2p + X(p). 

(3.42) 



Set 



and 



Z 2 {t) = j\puj 2 + 



2/x + Hp) 



)dx, 



V 2 (t) = j P{H 2 + L 2 )dx = J - [(jau Xi + F X2 ) 2 + {-puj X2 + F Xl ) 2 ]dx. 



Then it follows that for < r < i 



2" 



||V(*»|| 2(1 _ r) < C^(t)||p||L < C<p(t)(—)T& < C<p(t)r^, (3.43) 
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and 



||V«|| 2 + ||w|| 2 + ||divw|| 2 + ||(2// + A(p))5divit|| 2 
P{P? 



< C 



Z{t) + 



2/i + A(p) 



dx 



<C(Z(t) + l). 



(3.44) 



Now we estimate the four terms on the right hand side of (|3.42p . First, by the interpolation 
inequality and Lemma \2~2\ (I3.43P and (|3.44p . for < e < |, it holds that 



| I u 2 divudx\ < C7||divu|| 3 ||w||^ < C(Z(t) + 1)\\uj\\^ 2e \\Vu\\^ £) 

< C{Z(t) + l)Z(^^(i)^e( 1 -w-&) 

_ „ 2(l-2e) 2 l-£ 

< aip 2 (t) + C a Z(tf(Z(t) + 1)-L=3^£ W5T=3T 

< <V(t) + C a {Z{tf + l) 2+ T^ £ T^r^, 



(3.45) 



where and in the sequel a > is a small positive constant to be determined and C a is a positive 
constant depending on a. 

Next, one has 



- / F z divu 



1 



1 



dx\ 



2p J + X(p) J 2/x + A(p). 

1 f p2 f g , ^fc) \ 2/j + A(p) + pA'(p) 

2 7 ^2/i + A(p) + 2 / . + A(p)y' (2 M + A(p))2 1 



(3.46) 



P{P) 



2/i + A(p) 



dx < C 1 + 



\FY 



2/i + A(p) 



dx I , 



and 



Fdivu 



<ix| 



2/x + A(p) ' 2/x + A(p) 
1 t ( F P(p) \ P(p)(2p + X(p)) + pX'(p)P(p)-pP'(p)(2p, + X(p)) 

2.1 \2[i + \(p) 2p + X{p)J {2fi + X(p)) 2 



- / 'V2(. + A(p) 2,.+ A(rty w " V J 2l* + MP) J 
On the other hand, it holds that 

| - J 2F(ui X2 u 2xi ~ ui xi u 2x2 )dx\ < C j |F||Vu| 2 (ix. 



dx\ 

(3.47) 
(3.48) 



Substituting (l3T431) - (IQ8D into (l3Ti2|) yields that 
1 d 



2dt 



Z 2 (t)+^{t) 2 < a V (t) 2 +C a (Z(t) 2 +l) 2+ T^ £ —+C 



1 + 



2/i + A(/o) 



ds+ / \F\\Vu\ 2 dx 



Now it remains to estimate the terms 



(3.49) 



2/i + A(/o) 



dx and y |.F||Vu| 2 dx on the right hand side 



of (pT4"9|) . By Lemma E31 for e £ [0, \] and 77 = e, it holds that 

|[i^||2^ < Cni.Pll! +ml||V J P|| 1 i^[| J P||I (1 _ £) 



(3.50) 
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where s 



and the positive constant C is independent of m and e. 



m — e(l — e) 

Choose the positive constant e = 2~ m with m > 2 being integer in the inequalities (13.491) 
and (|3.50p . By the density estimate (|3.24p in Lemma 13.41 one can get 



ll*1h 



and 



ll^ll 



2(l-e) 



< 



(2^ + X(p)y- 2 \F\(2f, + X(p))2dx 

(2/i + A(p))"( 1 - £ )|F| 2 ( 1 - £ )(2/i + \(p)) 1 ~ £ dx 
\F\ 2 



(3.51) 



2(l-£) 



r. 



-Y 



V2^ + A( P ); 

< CZ(t) 



(2/x + A(p)) <= cfa: 



2(l-e) 



2 



i < cz(ty 



+ 1 



(3.52) 



sf3 



< cz(ty e~j^ + 1 = cz(ty 2^ + 1 < cz(t) s , 



ms/3 



where in the last inequality one has used the fact that ms = ;^r~77i~?i — )• 1 as m — )■ +00. 



Substituting (|3.43p with r 



to— e(1 — e) 

(133Tj) and (13321) into (133011 yields that 



||F|| 2m < C\z{t) + wh\\VF\\ 1 ^Z{t) 

< c 

Thus it follows that 



< C 



Z (t) + m^^jFi^) 1 "^^) 5 



Z{t)+m l *{™)T^y{t) l - s Z(ty 



(3.53) 



< 



2/i + A(p) ' 7 ( 2fi + x( p )) 1 -^T) 2(i + X(p) 

V2/i+A( P ); 1 1 

2m-3 

|F| 2 \ 2(m-l) 

2/i + A(p) 



)2(m-l) \F\ 1+ ^dx 



<-{S 



dx 



\F\ 2m dx 



2(m-l) 



< Z^^IIFH- 1 < CZ(t)^ [z(i) +mh{™)fe<p{t r -'Z{t) 

m (l-s)m (l-s)m (2 + s)m-3 



(3.54) 



< C 



Z(t) 6 + m 2 ^- 1 ' (— ) (^(t) m-i Z(t) ™- 

£ m 777, 2(l-a)m 2((2+s)m-3) 

Z(t) + m m ( 1+s )" 2 ( ) (/3-l)(m(l+s)-2) ^(t) m(l+<0-2 



< a<^(i) 2 + C, 

< a<^(i) 2 + C, 

where in the last inequality one has used the fact that ms 
m — > +00. 

Furthermore, it holds that 



Ul 2 _ o 1 1 — ms 

(1 + Z(t) 2 ) 2 + m{-)—^(l + Z{tff + ^+°)-i 

m(l-e) 2 _ 
m— e(l— e) 



1 with e = 2~ m as 



\F\\Vu\ 2 dx < ||Fl| 2m l|Vttl| 2 4 m < C\\F\\ 2m (\\dwu\\\r^_ + IM| 2 4 m ) 

2m-l V 2m — 1 2m-l/ 

<C||F|| 2m (|' 



2/i + X(p) 3^=1 



2 4m + 1 1 W 1 1 2 4m + 1 



(3.55) 
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Note that 



2/z + X(p) 



(2^ + X(p))—- 




dx 



\F\ C2m-l)(m-l) fa 



\F\< 



2(m-l) 



-fix 



2/i + A(p) 

and from ^ wdx = 0, Lemma 12.21 and (I3.43h . one has 



4(m-l) 

(2// + A(p)) 2m-3 / 
2m — 3 

|_P|2 \ 2(m-l) _J_ . 

dx] <C7||F||^Z(t)- 



(3.56) 



M|-4m <C|H|2 m(1 " 2E) HVwII^ 1 ,-;;' <CZ(tf"W e^(t) 

2m — 1 ' ' L 

m(l — e) Q 1 — £ 1 — £ Q 1 — £ 1 — £ 

< C2(' 3 - 1 ) m ( 1 - 2e )Z(t) m(l-aO <p(t) m(l-2e) < C Z (t) m(l-2e) y,^) m(l-2e) 



m(l-2e) 



Now substituting fl333j) . (l336jl and (l337|l into (13351 gives that 

o T 1 771 \—s -, I r o 1 — £ 1— £ 

|F||Vu| 2 <ix < C|Z(t) + m2(— )^( / 9(t) 1 ~ s Z(t) s J |l + Z(t) ^(t) 



1+: 



z(ty 



< c 



2 
1 + s 



(3.57) 



'z{t) + mH™)fe<p(t) x - a ztty 

Z{t) + Z(t) 3 + Z(tf~^^ <p(t)^™ + m^-)^^) 1 " 8 ^) 8 

+ m2( — )/3-l^(t) + ™(1-2 E ) Z(t) + m(l-2 e ) + m3 ( — )(/3-l)(m-l)(^(t) m-1 Z (t) rn-1 

< aipitf + C a [(1 + Z 2 {t)) 2 + (wk (™)Hz(ty 



1 777, 1— s o_i_ t . 1 — g 

+ ( m 2(^)FTZ(t) 2+s "^I 



m 



ms-1 
1 



^277 _|_ ( m 2( — )(/3-l)(m-l)^( i ) — 



2(m-l) 
m(s + l)-2 



< a^(t) 2 + C7 Q [(1 + Z 2 (t)) 2 + m(-) — (1 + Z{tf) 

m 2 „ r, j l-ms+(2ms-l)e yyj 2 „ o i 1-ms 

+ m( — )' 3 " 1 (1 + ^(t) 2 ) Z+ (l+ S )m(l-2 £ )-l+ e +m (_)/3-l(l + 2'(t) i )^ + '»( fl + i)-2 

Substituting (]3.54p and (I3.58P into (|3.49j) and choosing a sufficiently small yield that 



(3.58) 



~{Z 2 (t)) + \y{tf < C(Z(t) 2 + 1) 2+ i^ £ t^ + c[(l + Z 2 (t)) 2 + m (-)^T(l + Z(t) 2 ) 

Ul 2 „ r, j l-ms + (2ms-l)g j 2 „ o , 1-ms 

+ m( — )^"l (1 + Z(ty) Z+ (l+s)m(l-2e)-l+, + m a ( — )/3-l (1 + Z(i) 2 ) + m( s + l)-2 

(3.59) 

Note that lim m _j. +00 



large, one has 



[2 m (l — ms)] = 2, and so 1 — ms ~ 2e as m — > +00. Thus for m sufficiently 
1 — ms 2e 2e 



m(l + s) 



1 - 2e + m - 2 m - 1 - 2e 



and 



1 — ms + (2ms — l)e 



(1 - ms)(l -2e) + e 



3e - 4e 2 



(l + s)m(l-2e)-l+e (1 + s)m(l - 2e) - 1 +e (m + 1 - 2e)(l - 2e) - 1 + e 



< 4e. 
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Then f|3.59j) yields the following inequality for suitably large m, 



\jt {Z2{t)) + < Cm(^) A(l + Z(t) 2 ) 2 ^. 



Note that 



Z 2 (t) = j (piv 2 + 



-)dx 



P 2 (P) 



2p + X(p) 
<C |[^ 2 + (2^ + A( /0 ))(divn) 2 + 2 

< C(<f>(t) + J P 2 (p)dx) £ L 1 (0, T). 



)]dx 



Applying the Gronwall's inequality to (|3.60p and using (|3.61|) show that 

1 1 ,m _2_ 

(TT^WF " (i + ^(o))^ + Cm£{ 7 ] ^ ~ °" 

Then we have the inequality 



> 



provided that 



(1 + Z 2 (t))^ " 2(1 + Z 2 (0)) 4£: 



_ .m N _?_ 1 
Cm£( 7 )g -^ 2(1 + ^(0))^ - 



This condition, (|3.64|) . is satisfied if 



Cm 1+ Fi2- m(1 -Fi) < I, 
~ 2' 



since 



Z 2 (0) 



5\2 



TO 



< c 



l n oll# 2 (T 2 ) + 1 1 >°0 1 1 ^-3 1 1 t^o 1 1 1^ 2 (TT 2 ) + IIPolli/3( T 2 



< c. 



(3.60) 



(3.61) 



(3.62) 

(3.63) 
(3.64) 

(3.65) 



Now if j3 > 3, that is, 1 — > 0, then we can choose sufficiently large m > 2 to guarantee 
the condition (|3.65p . Consequently, the inequality (I3.63P is satisfied with (3 > 3 and sufficiently 
large m > 2. Then 

Z 2 {t) < 2 2 '™~ 1 (1 + Z 2 (0)) - 1 < C, (3.66) 



and 



[ tp{t)dt 
Jo 



< C. 



Thus the proof of Lemma 13.51 is completed. 

Step 5: Second order derivative estimates for the velocity: 



(3.67) 
□ 



Lemma 3.6 There exists a positive constant C independent of 5, such that 

sup f p(H 2 + L 2 )dx + f [ p(H Xl - L X2 ) 2 + (2/i + \(p))(H X2 + L Xl ) 2 dxdt < C. (3.68) 
te[o,T] J Jo J 
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Proof: Multiplying the equations, (|2.4p 1 and (|2,4p 2 . by H and L, respectively, summing the 
resulted equations together, and integrating with respect to x over T 2 lead to 

~ J p(H 2 + L 2 )dx + J fi(H xl - L X2 f + (2// + \{p)){H X2 + L Xl ) 2 dxdx 
= J p(H 2 + L 2 )divudx — J pudivu(L X2 — H Xl )dx 

- / ^2^ + A( P) ) ^(^^y)' + (^^) ) , ] di - + ( 3 - 69 ) 

— J \H(u X2 ■ VF + pu Xl ■ Vuj) + L(u Xl ■ VF — pu X2 ■ Vw)] dx 
+ J(2p + X(p))[(ui Xl ) 2 + 2ui X2 u 2xi + (u 2x2 ) 2 ](H X2 + L Xl )dx. 

Set 



and 



Note that 



Y(t) = (J p{H 2 + L 2 )dxj , (3.70) 
m = (J p(H Xl - L X2 ) 2 + (2/i + A(p))(tf, 2 + L,J 2 ^ 2 . (3.71) 

1 (|Vtf | 2 + \VL\ 2 )dx = J (Hi + H 2 X2 + L 2 X1 + L 2 X2 )dx 

= [ [(H Xl -L X2 ) 2 + (H X2 +L Xl ) 2 ]dx<-^ 2 (t). 
J A 1 

Thus it holds that 

\\V(H,L)\\ 2 (t) < Ci/>(t), Vt G [0,T]. (3.72) 
Then it follows from the elliptic system 

puj Xl + F X2 = pH, -poj X2 + F X1 = pL, 

that 

\\V(F,u)\\ p <C\\p(H,L)\\ p , VKp<+oo. (3.73) 

Furthermore, since ^ (p>w Xl + F X2 )dx = 0, by the mean value theorem, there exists a point 

x* £ T 2 , such that (poo Xl + F X2 )(x*,t) = 0, and so H(x*,t) = 0. Similarly, there exists a point 
x*, such that L(x' Jf ,t) = 0. Therefore, by the Poincare inequality, it holds that 

\\(H,L)\\ P <C\\X7(H,L)\\ 2 , Vl<p<+oo, (3.74) 

where C may depend on p. 

Now we estimate the right hand side of (|3.69|) term by term. First, by the Holder inequality, 
(|3.74p and the density estimate (|3.24p . it holds that 



J p{H 2 + L 2 )divudx| = | J p(H 2 + L 



2P + \{P) (o 7 r\ 

< \\^p(H,L)\\ 2 \\(H,L)\\ 4 \\ ^ P[ Zu * ^ CY(t)ma + \\F\k). 

2p + A(p) 
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Note that 
\\(F,u)\U<C(\\V(F,u)h + \\(F,u) 



< C 



/ r' 2 f 2 \ 1 / r \ l i r i ( 3 - 76 ) 

||V(F,o;)||| + ( J 2 x{ ) dx) 2 ( J (2p + X(p))dxj 2 + || W || 2 J < C[Y(t) + 1 



where in the last inequality one has used the estimate (|3.43j) with r = \ and the estimate (|3.66j) . 
Substituting (|3.76|) into (|3.75|) yields that 

| J p{H 2 + L 2 )divudx\ < CY(t)i/t{t)(Y(t) + 1) < ai) 2 {t) + C a (Y(t) + l) 4 . (3.77) 

Second, direct estimates give 



| — J poodivu(L X2 — H Xl )dx\ < p (^j (L X2 — H Xl ) 2 dx^j (^j oj 2 (divu) 2 dx 

< a^(t) + C a J u 2 {d\vufdx < a ^t) + CJu\\l\\ ^±^ \\l ^ 

< a^(t) + C a \\u\\l(l + \\F\\l) < m/> 2 (t) + C a (Y(t) + l) 4 . 



Similarly, one has 

| - / ,(2, + AW) [Hjj^Y + + L n )M 

< a j (2/1 + \{p))(H„ + L n fdx 

< a^ 2 (t) + C a (l + \\F\\i) < «V 2 W + C a (Y(t) + l) 4 . 



Next, 



| - J [H(u X2 ■ VF + fiu Xl ■ Vu) + L(u Xl ■ VF - pu X2 • Vw)] dx\ 

< cj \{H,L)\\Vu\\V{F,uj)\dx (3.80) 

< C||(iJ,L)|| 8 ||V«|| 2 ||V(F,a;)||| < C\\V(H, L)\\ 2 \\p(H, L)\\s, 



where one has used the fact that 



||Vn|| 2 < C(||divw|| 2 + 1Mb) < C(|| f + P ,[ P \ h + < C. 

Ip + A{p) 



Note that 



\\p(H,L)\\s = U p?\(H,L)\*dx\ = U y/p\(H,L)\\(H,L)\*p-rdx 
< \\^(H,L)\\l\\(H,L)\\l\\p\\i < CY(t)l\\V(H,L)\\l 
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(3.81) 
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It follows from (01801) and (ET8D that 

| — J [H(u X2 ■ VF + pu xi ■ Vw) + L(u Xl ■ VF — pu X2 ■ Vw)] dx\ 



(3.82) 

<CY{t)i\\V{H,L)\\J <CY{t)i^{t)f < a^(t) 2 + C a Y{tf. 



Moreover, 



| J{2p + X(p))[(ui xi ) 2 + 2u lx2 u 2xi + (u 2x2 ) 2 ](H X2 + L Xl )dx\ 

< a^(t) 2 + C a J(2fx + X(p))[(u lxi ) 2 + 2u lx2 u 2xi + (u 2x2 ) 2 ] 2 dx 

<am 2 +C a \\2p + \(p)\\ 2 \\Vu\\i ( 3 - 83 ) 

< cn/^) 2 + C a (||divu||t + ||w||f) < am 2 + C a {\\{F,u>)\\j + 1) 
<ai;(tf +C a (\\V(F,u)\\i +1) <am 2 + C a (l + Y(t)) 4 . 

5 

Substituting the estimates (|3.77|) - (|3.79p . (|3.82p and (|3.83|) into (|3.69p . one can arrive at 

\jt {Y2{t)) + ^ {t) ~ 5a ^ {t) + Ca{1 + y2(t))2 - (3 - 84) 

Choosing 5a = ^, noting that Y 2 (t) = (f 2 (t) G L 1 (0,T), and then using Gronwall's inequality 
yield that 

Y 2 (t)+ [ iP 2 (t)dt<Y 2 (0) + C. (3.85) 
Jo 

Now we calculate the initial values y 2 (0). By the approximate compatibility condition (13. 21) . 
one has 

£ p su S - VP 5 = ^ Q g, with <? G L 2 (T 2 ). 
On the other hand, it holds that 

C pS u 5 = pAu s + V((p + A(^))div«g) = pAu s + V(F 5 - /xdivug + P 5 ) g 
= [/iV(div^) -/iVx(Vx «*)] + V(F 5 - pdivu 5 + P 5 ) 

where Fq = (2p + X(pq))(Hvuq — Pq and similarly one can define Wq, Lq,Hq, Vx denotes the 
3-dimensional curl operator, and 

v x iy x 4) = (d X2 J Q ,-d Xl 4,o) 

is regarded as the 2-dimensional vector (8 X2 ujq, — c^Wg)*. 
Thus 

£,gug - VP* = VF* - ^wg, -d^)* 

= (i^ - pd X2 u; 5 , F$ X2 + pd Xl 4f = p 5 (L s , H^f. 



(3.87) 



Therefore 

^g = p s (L 5 ,H s ) t . (3.88) 

Consequently, it holds that 



3^(0) = \Wp 5 o(L s Q ,H s )\\l = \\ZfLgWi < C. (3.89) 
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This, together with (|3.85j) . shows that 

Y 2 (t) + f ^ 2 (t)dt < C. (3.90) 
Jo 

This completes the proof of Lemma 13.61 □ 

Remark 3.1 Similar to the derivation of (j3.87|) . one can get that for any t £ [0, T\, 

C p u-VP(p) = p(L,H) t . 
Then it follows from the momentum equation (|l.ip 2 that 

u t = (L,HY -u-Vu. (3.91) 
The above identity can also be obtained directly from (|2.1|) . 

Step 6. Upper bound of the density: We are now ready to derive the upper bound for the 
density in the super- norm independent of 5, which is crucial for the proof of Theorem 11.11 as in 
~j3l. First, we have 



Lemma 3.7 It holds that 

rT 



im^WldtKC. (3.92) 



Proof: By (|3.73j) with p = 3, one has 



T \\V(F,u;)\\ldt <C \\p(H,L)\\ 3 3 dt = C J J p 3 \(H, L)\ 3 dxdt 
= cjj ^\(H,L)\\(H,L)\ 2 pUxdt 

<C j ||^^)||2||(^,i)||i||p||fo^ 

<C [ \\V(H,L)\\ 2 2 dt <C [ ^ 2 {t)<C, 
Jo Jo 

which, combined with the estimates in Lemma 12.31 yields that 



(3.93) 



[ T IK^o;)!! 3 ^^ < f T 11(^)11^,3^ < C. (3.94) 
Jo Jo 

The proof of Lemma 13.71 is finished. □ 
With Lemma 13.71 in hand, we can obtain the uniform upper bound for the density. 

Lemma 3.8 It holds that 

p(t,x) < C, V(t,x) G [0,T] x T 2 . (3.95) 
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Proof: From the continuity equation (jl.ip -p we have 

6(p) t +u-Vd(p) + P{p) + F = 0, (3.96) 

where 6(p) is defined in ()3,22l) . 

Along the particle path X(r;t,x) through the point (t,x) £ [0, T] x T 2 defined by 

dX(r;t,x) , ,-», .. 

V =u(T,X(T:t,x)), 

dr v v (3.97) 

^ X(r;t,x)\ T=t = x, 
there holds the following ODE 

^e(p)(r,X(r;t,x)) = -P(p)(r,X{r;t,x)) - F(r, X(r;t,x)), (3.98) 
dr 

which is integrated over [0, t] to yield that 

0(p)(fi, x) - 6(po)(X ) = - f\p(p) + F)(t, X(t; t, x))dr, (3.99) 

J o 

with Xq = X(t; t, x)| r= o. 

It follows from ([3^99]) that 

o(t x) 1 /"* 1 ^ /"* - 

2fihx^{+-^(t,x)+ P(p)(T,X(r;t,x))dr = -p (X f- F(r,X(r;t,x))dr. (3.100) 
po(^o) P Jo P Jo 

So 



Po(X ) P 
which implies that 

p(t,x 



2pln-^ z <l\\p \t+ /V(^-)||oodr<C, (3.101) 

Jo 



< C. 



Po(X ) 
Therefore, we have 

p(t,x) < C, V(t,x) £ [0,T] x T 2 . (3.102) 
Hence the Lemma is proved. □ 
As an immediate consequence of the upper bound of the density, one has 

Lemma 3.9 It holds that for any 1 < p < oo, 

/ (HdiHlL + \\V(F,u)\\l)dt < C. (3.103) 
J o 



Proof: First, note that 

fT 



[ \\divu\\ldt<C [ m\lo + \\P(j>)\\lo)<tt<C. (3.104) 

J •/ 



Then for any 1 < p < oo, 

T \\V(F,u)\\ldt <C I \\p(H,L)f p dt 



o 



(3.105) 



<C [ \\(H,L)\\ 2 p dt<C [ \\V(H,L)\\ 2 2 dt<C. 
Jo Jo 

Thus Lemma 13.91 is proved. □ 
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4 Higher order estimates 



With the approximate solutions and basic estimates at hand, we can derive some uniform esti- 
mates on their higher order derivatives easily as in \25\ [22| [23] . We start with estimates on first 
order derivatives. 



Lemma 4.1 It holds that for any 1 < p < +oo, 



sup ||(Vp,VP(p))(V)|| P + 
te[o,T] Jo 



\Vu\liodt < C. 



(4.1) 



(4.2) 



Proof: Applying the operator V to the continuity equation (JTTTJ) 1 , one has 
(Vp) t + Vu-Vp + u- V(V» + Vpdivu + pV(divu) = 0. 
Multiplying the equation (|4.2|) by p|V p\ p ~ 2 V p with p > 2 implies that 

{\Vp\ p ) t + dw(u\Vp\ p ) + (p-l)\Vp\ p divu + p\Vp\ p ~ 2 Vp-(Vu-Vp)+pp\Vp\ p ~ 2 Vp-V{dwu) = 0. 

Integrating over T 2 gives 
d 



(4.3) 



dt 



l|Vp||£ 



-0-1) J \Vu\ p divudx - p J \Vp\ p ~ 2 Vp- {Vu-Vp)dx - p J p\Vp\ p ~ 2 Vp ■ V(divu)dx 
< (p- l)||divu|| 0O ||V/>||£ + PllVtilloollVpllP+pllplUllVpll^UVdiHIp. 



(4.4) 



This implies that 



^llVpUp < C7[||V«||oo||Vp|| p + ||Vdivu|| p 



< C 

< c 



|V«|I-I|V* + ||V(£±M)||, 
VnHoo + Halloo + l) || Vp|| p + ||VF|| P 



(4.5) 



By Remark 13 -1| one has 



C p u = VP(p)+p(L,H) t . 



Thus the elliptic estimates and (|3.74p yields that for any 1 < p < oo, 

\\V 2 u\\ p <C[\\VP(p)\\ p + \\p(L,H)\\ p ] 

<C[\\Vp\\ p + \\(L,H)\\ p ) <C[\\Vp\\ p +\\V(L,H)\\ 2 } 

By Beal-Kato-Majda type inequality (see [23]- [25] or [51]), it holds that 



(4.6) 



(4.7) 



|Vlt||oo < C(||divit||oo + IMloo 
< Cdldivulloo + \\cj\\ 

oo 



ln(e+ ||V 2 n|| 3 ) 

ln(e + ||Vp|| 3 ) + C(||divu|U + IMU) ln(e + ||V(iJ, L)\\ 2 ) 



(4.f 
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The combination of (|4.5p with p = 3 and (j4.8j) yields that 



d 



^||Vp|| 3 < C[(||divu||oo + H^Hoo) ln(e + ||V(F,L)|| 2 ) + ||F||oo + lj ||Vp|| 3 
+ COldivulloo + IMloo) l|Vp|| 3 ln(e + ||Vp|| 3 ) + C||VF|| 3 . 

By the estimates (|3.94p . (|3.104|) . (|3. 105|) and the Gronwall's inequality, it holds that 

sup ||Vp|| 3 <C, 
te[o,T] 



which, together with (J3J3D, (|3.1U4I) . (f4T7|> and (|4Tg]l . yields that 

HVulliLcft < c 



o 



(4.9) 



(4.10) 



(4.11) 



Therefore, by (|4.1ip . Lemma [3.71 Lemma [3. 9 1 and Gronwall inequality, one can derive from fj4.5j) 
that 

sup ||Vp|| p < CdlVpollp + 1), Vp€[l,+oo). (4.12) 

te[o,T] 

Thus the proof of Lemma |4. II is completed. □ 



Lemma 4.2 It holds that for any 1 < p < +oo, 



sup 

te[o,T] 



|«(t,-)||oo + ||V^|| p +||(p t ,P t )|| p +||(p t ,P(p) t )|| if i + ||(p,«)|| H2 +/ |H|^ 3 dt<C. (4.13) 

•/ o 



T 



Proof: By L 2 — estimates to the elliptic system (|4.6p . one has 



sup |k||^ 2 < C sup (\\VP(p)\\ 2 + \\p(H,L)\\ 2 ) 
te[o,T] te[o,T] 

<C sup (||VP(p)|| 2 + ||v^(#,L)|| 2 ) <C. 

te[0,T] 



(4.14) 



It follows from the Sobolev embedding theorem that 



sup \u(t,x)\<C, sup ||V«|| P < C, VI < p < +oo. 

[0,T]xT 2 te[0,T] 



(4.15) 



Due to (jl.ip i . one can get pt = —u-Vp — p divu and Pt = —u-S/P — pP'(p) divu, which, together 
with the uniform upper bound of the density and the estimates in Lemma l4.1l and (|4.15p . yields 
that 

sup \\(pt,Pt)\\p<C, Vpe[l,+oo). (4.16) 
te[o,T] 

Applying V 2 to the continuity equation (j 1 . 1 [> x . then multiplying the resulted equation by V 2 p, 
and then integrating over the torus T 2 , one can get that 



d_ 
di 



|V 2 p|| 2 < C 

< c 



|Vi*[|oo|[v 2 p[|l + ||Vp|| 4 ||v 2 p||2||v 2 «||4 + l|p||ool|v 2 p|| 2 ||v 3 u|| 2 

;i|Vu||oo + l)||V 2 p||l + ||V 3 u||| + l 



(4.17) 
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Similarly, 



j t \\V 2 P(p)\\ 2 2 <C 



iVwIloo + l)||V 2 J P(p)||§ + II V 3 n||? + 1 



(4.18) 



Note that (|4.6|) implies that 

£ P (V«) = V 2 P(p) + V[p{H, L)] + V(VA(p)divn) := $. 
Then the standard elliptic estimates give that 

|M|ff3 < C 



< C 



U\\ H 1 + ||$||2 

u \\ m + ||v 2 P( P )|| 2 + ||p|U||v(#,L)|| 2 + ||Vp|| 4 ||(ir,L)||4 

+ ||V 2 p||2||divtt||oo + ||Vp||4||V 2 «|| 4 l , 



(4.19) 



and 



Consequently, 



|V 2 n|| 4 < C||VP(p)|| 4 + \\ P (L,H)\U < C(l + ||V(P,L)|| 2 ). 



u\\h 3 < C 



1 + ||V 2 P(p)|| 2 + ||V(P,L)|| 2 + ||V 2 p|| 2 ||divu| 



Substituting (ETT9]) into (j4"T71) and (j4~T8|) yields that 

|||(V 2 p, V 2 P(p))|| 2 < c[(||Vu||L + l)||(V 2 p, V 2 P(p))|| 2 + ||V(P,L)|| 2 + 1 



(4.20) 



(4.21) 



Then the Gronwall's inequality yields that 



C 



T 



||(V 2 p,V 2 P(p))|| 2 (t) < (J|(V 2 po,V 2 P )|| 2 + C J (\\V(H,L)U + l)dt)e Jo 

<c, 

which also implies that 

sup (\\(j>,P(j>))\\H* + \\(jH,P(p)t)\\m)+ f T \\u\\%sdt<C. 

tG[0,T] J 

The proof of Lemma 14.21 is completed. 



Lemma 4.3 It holds that 



fT 

sup || v / Mlll(*)+ / Klllri^ <C 

!GfO,Tl JO 



te[o,T] 

Proof: The momentum equation (ILlh o can be written as 

pu t + pu • Vu + VP(p) = £ p it := pAu + V((p + A(p))divu) 
Applying c?t to the above equation gives that 



(4.22) 

(4.23) 
□ 

(4.24) 
(4.25) 



pu tt + pu ■ Vu t + VP(p)t = pAu t + V((p + A(p))divu 4 ) - p(U t - p t u ■ Vu - pu t ■ Vu + V(A(p)tdivu) . 

(4.26) 
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Multiplying the equation (|4.26p by ut and integrating the resulting equation with respect to x 
over T 2 imply that 

\dt S p ^ 2dx+ f (^i v ^i 2 +(^+ A (p))i div ^i 2 )^ 

= - J VP(p)( • u t dx - j Pt \u t \ 2 dx - j Pt (u ■ V«) • u t dx (4.27) 



p(ut ■ Vu) ■ utdx + I V(\(p)tdivu) ■ utdx. 



Notice that 



- j VP{p) t • u t dx = J P{p) t d\\u t dx 

<^J \divu t \ 2 dx + C J \P t \ 2 dx <^J \divu t \ 2 dx + C, 

J pt\u t \ 2 dx = j div(/9n)|n t | 2 (ix = — 2 j p(u ■ Vu<) • u t dx 

<^J iV^l^x + CIHlLllv^llLll^ll^^ / \Vu t \ 2 dx + C\\^pu t \\l 



and 

V (X(p)tdivu) • utdx = — I \(p)tdivudivutdx 



\dwut\ 2 dx + C\\X(p)t\\l\\dwu\\l < | J \dwu t \ 2 dx + C. 



(4.28) 



(4.29) 



(4.30) 



p t {u ■ Vu) • utdx = J div(pu)[(u ■ V«) • ut]dx = — J pu -V[(u ■ Vu) ■ ut]dx 

< llp||oo|k||Ll|Vnt|| 2 ||Vn|| 2 + hlloollVplloollv^^lbdlVnlH + ||«||oo||V 2 «|| 2 ) 
<| J \Vu t \ 2 dx + C{\\^pu t \\ 2 + \\(Vu,V 2 u)\\ 2 + \\Vu\\i) 

<%J \Vu t \ 2 dx + C{\\^u t \\ 2 + l), 

p(u t • Vu) • u t dx\ < WVuWooW^utWl (4.31) 



(4.32) 



Substituting the above estimates into (|4,27p and then integrating with respect to t over [0, t] 
yield that 

fwVutgdtK \\J44m\l + C [ (||V«||oo + l)||^ll^ + C. (4.33) 
Jo Jo 

By the compatibility condition (|3.2p . it holds that 



l 0i 

•4 



thus we have 



p s uf(0)\\ 2 < \\^kg\\ 2 2 + ll^llooll^llLllV^H 2 < C, 
Po 



which, together with (|4.33p and the Gronwall's inequality, yields that 



sup \\Jpu t \\l{t) + / ||Vu t |||di < C. 
te[o,T] Jo 



T 



By dMO, for any 1 < p < +00, 

fT r-T 



[ \\ut\\ 2 p dt< [ (\\(H,L)\\ 2 p + \\u\\l\\Vu\\ 2 p )dt 
Jo Jo T 

< f (\\V(H,L)g + \\u\\l\\Vu\\l)dt<C. 
Jo 



Therefore, one can arrive at 



f \\utf H idt < a 

Jo 

Thus the proof of Lemma 14.31 is completed. 



Lemma 4.4 It holds that 

sup \\(p t ,P(p) t ,X(p) t )\\ H i(t)+ f \\{ptt,P(p)ttA{p)tt)\\ldt<C. 
te[o,T] Jo 

Proof: From the continuity equation, it holds that p± = —u ■ Vp — pdivu and pu 
u ■ Vpt — ptdivu — pdivut, and thus 



te[o,T] 



sup ||Vpt||2(i) < sup ||Vp||4||Vu||4 + || ||V 2 p|| 2 + ||p||oo||V 2 u| 



te[o,T] 



< C. 



and 



T r T 

WpttWldt < 
Jo 



il|Vp|| 2 + IMILIIVaH 2 . + ||pt|| 2 ||Vn|| 2 + llpllSollVutlH 



(It 



<C [ T (\\u t f H1 +l)dt<C. 
Jo 



Similarly, we have 

sup ||V(P(p) t ,A(p) t )|| 2 (t)+ [ T \\(P(p)u,\(ph)\\ 2 2 dt<C. 
te[o,T] Jo 

Thus the proof of Lemma 14.41 is completed. 
Lemma 4.5 It holds that 

H1 +t\\ u \\i 3 +t\\(p tt ,p(p)uA(ph)\\i + \\(p,p(pmw^ 



sup 

te[o,T] 



T 

"/ ' 
10 



H 2 (t) + \\ u \\h 4 



dt < C. 
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Proof: Now multiplying the equation (|4.26p by uu and then integrating with respect to x over 
T 2 yield that 



IIVP"tt|li(*) + ^ J (m|V^| 2 + (/i + \{p))\d\\u t \ 2 ^dx = i y A(p) t |divitt| 2 <£c 

- / (VP, + ^ + • V. + P u ■ Vn, + p., • Vn) • u tt dx + / V(A(p) t divn) • n t ^*. 



(4.40) 



Note that 



V(A(p)tdivu) • uudx 



d 



X(p)tdivudivuudx 



— j \(p)tdivudivutdx + ^ (A(p)t|divu| + X(p)ttdivudivu t ^dx. 



Substituting the above identity into (|4.40j) yields that 

\\Vpu t t\\l(t) + ~ J (p\Vu t \ 2 + {p + A(p))|divn t | 2 + \(p) t divudivu^jdx = ^ J X(p) t \divu t \ 2 dx 
— J (VPt + PtUt + Pt.u ■ Vn + pu ■ Vn t + pu t ■ Vn) • u tt dx + J X(p) tt divudivu t dx. 



Note that X(p) satisfies the transport equation X(p)t 
holds that 



J A(p)t|divnt| 2 cte| 



n • VA(/3)|divnt| 2 (ix 



(4.41) 

-n • VA(p) — pA'(p)divn, and then it 
pX' (p) divn | divnt 1 2 dx \ 



X(p)divutu ■ V (divu t )dx + - J (A(p) — pA'( / o))divn|divn t | 2 (ix| 
< C||A(p)n||oo||divnt||2||V(divnt)||2 + C\\X(p) - / oA / ( / 9)|| 00 ||Vn||oo||divnt| 



< C 1 1 div7/f 1 1 2 1 1 ^ (divitf ) 1 1 2 + CUVnlloolldivn 



* Ha- 



lt; follows from (p~26j) that 

C p u t = pu tt + PtUt + ( P u ■ Vu) t + VP(p) t + V(A(p) t divn). 
Then the standard elliptic estimates show that 



(4.42) 



|V 2 n t || 2 < C 



llv^lUIIVMilk + llptlWKIk + \\pt\U\ 



|Vn|| 4 + ||p||oo||«i||4||Vn||4 



+ ||pw||oo||Vnt|| 2 + \\VP{p)th + ||VA(p) t ||2||divn|| 00 + ||A(p) 4 || 4 ||V 2 n||4 



< C 

< c 



\\\fpu t t\\2 + \Wt\U + 1 + l|Vnf|| 2 + ||divn||oo + ||V 2 n|| 4 
\\y/puuh + IKIk + 1 + ||Vn f || 2 + ||divn||oo + ||V 3 n|| 2 



(4.43) 



Substituting (031 into (Q2l) yields that 



A(p) t |divnt| 2 dx| < -|| ^pu tt \\ z 2 + CdlVnlU + l)\\Vu t \\ z 2 + C(\\u t \\j + ||V d n||^ + ||Vn||^). 

(4.44) 



At the same time, it holds that 

yP{p)t • u tt dx = j P(p) t d\\uttdx = J P(p) t d\\u t dx - J P(p)udivu t dx 
d f 

- ~dt\ P (p)t divu t dx + ll-PO^ttll! + ||div7Xt III, 



(4.45) 
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J ptUt ■ uttdx = J 



> t |2 d 



Kl 2 . 

pt—r-dx 



Ptt—^—dx, 



while 



M 2 , 
p tt ——dx 



div (pu) J Ut ^ dx 



(pu) t ■ Vu t ■ u t dx 



< WVpWocWVputhWutMVutiu + Hooll/otlUlklUIIViitllj 



< c 



HutlUllVtitlU+KlkllVtitll; 



< c 



t||V 2 u t || 2 + ||n 4 || 4 ||V^|b 



<§llv^ttll2 + C 
Moreover, it follows that 

ptu ■ Vu ■ uudx 
d 



< C\\u t \\ 4 W^/puuh + |MU + 1 + ||Vut|| 2 + Hdivulloo + \\V s u\\ 2 

l + l + \\Vut\\ 2 2 + \\d\vu\\l + \\VM 2 2 



1 



■— / ptu ■ Vu ■ utdx + / puu-Vu-utdx / ptUfVu-utdx 

d J f J 

- — / p t u- Vu- u t dx + ||ptt|| 2 |M|oo || V«||4 



ptu • Vut ■ utdx 



+ 



[||V«||4 + 



I^Hooll V^tlbll^tlU 



< 



d 

dt 
d 

"dt 



-~Jt I PtU ' Vu ' Utdx + c H'°* t ll 2 ll u *ll 4 + M2 + Il v ^ll2 



ptu ■ Vu ■ utdx + C 



W1 + KII2 + 1^111 



J pu ■ Vu t ■ u tt dx < Hv^ttlbllv^lloollVutl^ < gUVP^ttlll + C|l Vn tll 



pu t ■ Vu ■ uudx ||- v /7wtt||2|lv / Pll°°ll Vll ll< 



< o \\VPUtt\\ 2 2 + C\\u t \\ 



■1; 



and 



X(p)t t divudwu t dx < ||A(p) tt || 2 ||Vn|| 00 ||divu t ||2 < - ||A(p) tt ||| + || VuH^ ||divut \\l 
Collecting all the above estimates and substituting them into (|4.4ip yield that 
\\\^PUtt\\l{t) + jG{t) 



< C 



\\(ptt,P(ph,X(p)u)\\ 2 2 + 



+ ||V 3 u||| + (||V W ||^ + l)(||V^||^ + l) 



where 



(4.46) 



(4.47) 



(4.48) 



(4.49) 
(4.50) 

(4.51) 



(4.52) 



/\u P 
(^p\ Vu t | 2 + {p + A(p)) |divu t | 2 + X(p) t divudivu t - P(p) t divu t + Pt—i — I- Ptu • Vu ■ u^j dx. 

(4.53) 



Note that 



A(,9)tdivndivti t (ix| < \\\{p) t ||4||divti|| 4 ||divu t || 2 



< ^\\Vu t \\i + C\\X{p) t \\l\\dwu\\l< "-\\VutWi + a 



2 / Ml 
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J P(p) t dwu t dx\ < |||Vn t ||| + C\\P(p) t \\l < |||V«t||l + C, 



Ut j i 



U t 



I div (pu) - <ix| = | y pu ■ V ut ■ utdx\ 
< ll^|| 2 ||^||oo||Vn t || 2 < ^IIV^IH + C, 



and 



y ptu ■ Vu ■ u t dx\ = | j div(pu)(u ■ Vu ■ Ut)dx\ = \ J pu - V(tt ■ Vu ■ ut)dx\ 

< llv^llallv^lloodlVulU + |M|oo||V 2 u|| 2 ) + ||p|«| 2 ||oo||Vt*t||2||VM||2 

< #l[v«t||i + c. 



Therefore, it holds that 



Ci(||V« t |||-l)<G(i)<C(||V« t ||i + l), 



(4.54) 



for some positive constants C,C±. 
Now from (16.91). we can arrive at 



\\\V^u\\l{t) + j t G{t) 

< c\\\{p tu P{p) tu \{p)tt)\\l + + ||V 3 ni + (||Vn||L + l)(G(t) + 1) 



(4.55) 



Multiplying the above inequality by t and then integrating the resulting inequality with respect 
to t over the interval [r, t±] with r, t\ G [0, T] give that 



t\\^putt\\l(t)dt + ^ CtG(t) + c 

-tl 



dt 



+C 



\\(pu,p(p)u, Hphm + Kill + \\vM\i + G(t) 



(4.56) 



dt. 



It follows from Lemma [4.31 and (|4.54p that G(t) G L 1 (0, T). Thus, due to [6], there exists a 
subsequence t\. such that 

Tfc -> 0, T k G(r k ) ->• 0, as fe ->■ +oo. (4-57) 
Taking r = in (|4.56|) . then A; — >■ +oo and using the GronwalPs inequality, one gets that 

sup [t||Vut||!(i)] + f t\\^pu u \\l(t)dt < C. (4.58) 
te[o,T] Jo 



Note that (j4T4"3|) implies that 



sup [t||Ut,P(p) tt ,A(p) tt )||i(t)]+ / t\\V 2 u t \\ 2 2 (t)dt < C. 
te[o,T] Jo 



(4.59) 



It follows from (|33T]) that 



V = Vut- V(u- Vu). 
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Consequently, it holds that 



sup [t\\V(L,Hff 2 (t)] <C, 
te[o,T] 



which, together with ()3.74p . implies that 



sup [tW^HfWHt^^C sup [t\\V{L,H)%{f)]<C. 

te[0,T] t€[0,T] 



Therefore, it follows that 



sup [t\\u t \\l{t)] < C sup [t||(L,fl-)*||l(t)+t||«.Vu||i(t)] <C. 
te[o,T] te[o,T] 



So one can infer further that 



sup [t\\ut\\ 2 H i(t)] + / iK||| 2 (i)dt<C. 
te[o,T] Jo 

Applying d XjXk , j, k = 1, 2, to flUK gives 

(Px jXk )t + U ■ V(p XjXk ) + Ua. ja . fc • Vp + « Xj • Vp X)fe + U Xk ■ Vp Xj 

ivu) Xfe + p Xk idivu) x . + /?(divu) XjXj , 



(4.60) 



(4.61) 



(4.62) 



(4.63) 



0. 



Multiplying the above equation by q\ V 2 p\ q 2 p XjXk with q > 2 given in Theorem 1 1 . 1 1 and summing 
over j,k = 1,2 give that 



(iVVHi + divHVVl 9 ) + (<7 - l)|VVl"div U + q\V 2 p\ q - 2 p 



XjXk 



+u Xj ■ Vp Xk + u Xk ■ Vp Xj + p Xj (divu) Xk + p Xk (divu) Xj + p(divu) XjXk 
Integrating the above equality with respect to x over T 2 leads to that 



u XjXk ■ Vp 
0. 



d 



Thus one can get 



d l|V 2 p||, < C 



di 



< C 



||Vn|| 00 ||VV|| g + ||Vp|| 23 ||V 2 «|| 23 + 



llVulloollV 2 ^!, + ||V 2 u|| w i„ 



»l|V 3 u||, 



(4.64) 



where q > 2. Similarly, one can obtain 



d 



^||V 2 P|| g < ^[llVnlUllV 2 ^^ + \\V 2 u\ 



w 1 ^ 



(4.65) 



Apply d Xi with i = 1, 2 to the elliptic system C p u = put + pu ■ Vu + VP(p) to get 
£ p u Xi = -V(X(p) Xi divu) + p Xi u t + pu Xit + p Xi u ■ Vu + pu Xi -Vu + pu - Vu Xi + VP(p) x% := *5>. 
Then the standard elliptic regularity estimates imply that 



\\Vu\\ W 2, g < C 

< C 

< c 



l|Vu||, + ||*||,j 
1 + (||Vw||oo + l)ll(V 2 p, V 2 P)||, + ||Vu| 



W 1 -* 



(4.66) 



1 + (||Vx*||oo + l)||(V 2 p, V 2 P)|| 9 + \\u\\ H3 + \\u t \\ m + \\Vut\\ c 
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Thus it follows from iOl . (H~65l) and (IOTP that 

A||(V 2 p, V 2 P)|| g < C[l + (||Vu[U + 1)||(VV V 2 ^)!!, + ||n|| H 3 + II^Hhi + ||Vu t ||J . (4.67) 

Note that Lemma 12.21 implies that 

f' T f T l' T i 

/ \\Vu t \\ q (t)dt <C \\V 2 u t \\ 2 (t)dt < C sup [Vt\\V 2 uth(t)] / t'ldt < C. 
Jo Jo *e[o,T] j o 

Therefore, it follows from (|4.67p and the Gronwall's inequality that 
\\(V 2 p,V 2 P(p))Ut) < (\\(V 2 p ,V 2 P(p ))\\ q 



rt C I 

+C J (1 + \\u\\ H3 + ||ut|| H i + \\Vut\\ q )ds)e Jo 



t 

C / (HVulUW + lJds 

<C, 



(4.68) 



which then gives 

sup \\(p,P(p))\\w*«<P)<C. (4.69) 
te[o,T] 

So the proof of Lemma l4"31 is completed. □ 
Lemma 4.6 It ZioZck f/iai /or any < t < T, 

sup [i 2 ||v / P«tt||IW + * 2 ||«t||^ +* 2 II«II^3,,] + / t 2 \\Vu u \\ 2 2 {t)dt < C. (4.70) 
te[o,T] j o 

Proof: Applying to the equation (|4.26p gives that 



puttt + pu ■ Vu tt - CpU U = -Vptt - Ptt(u t + u ■ Vu) - 2pt(u u + u t -Vu + u- Vu t ) 
-2pu t ■ X7u t - pu tt ■ Vu + 2V((A(p)) t divu t ) + V(A(p) ti divu). 



(4.71) 



Multiplying the equation (|4.7ip by utt and integrating the resulting equation with respect to x 
over T 2 yield that 

\~dt J p \ Utt \ 2dx + J H Vu d 2 + (A* + \(p))(d\\uu) 2 dx = J pudivuttdx 

- / Pu{ut + « • Vu) • u tt dx - 2 / ptiutt + t* • V« + « ■ Vu t ) • u tt dx - 2 / pu t • Vu, • u tt d, 

- J pu tt ■ Vu ■ u tt dx - 2 / A ( p) t divn t divn t ^x - / A(p) ti div.divn^x. 

Multiply the above equality by t 2 to get that 

2c^(* 2 / P\ Utt \ 2dx ) ~ t [ pWtt\ 2 dx + t 2 J p\Vu tt \ 2 + {p + \{p)){divutt) 2 dx = t 2 J Pttdivuttdx 
-t 2 J p u (ut + u ■ Vu) ■ uttdx - 2t 2 J pt(u a + u t .Vu + u. Vut) ■ u tt dx - 2t 2 J pu t • Vu t • u tt dx 

7 

— t 2 / putt • V u • uttdx — 2t 2 / X(p)tdivutdivuudx — t 2 / \(p) t tdivudivuttdx := 
J J i=i 

(4.72) 
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Clearly, 



|/i| < at 2 \\divu tt \\ 2 2 + C a t 2 \\P tt \\l 



Now we estimate I 2 , which is a little more delicate due to the absence of estimates for u tt - First, 
rewrite I 2 as 



h=t 2 J dw(pu) t (L, Hf ■ u tt dx = -t 2 J (jm) t ■ V [(L, Hf ■ u tt ) dx 

= -t 2 J put-V [{L, Hf ■ u tt ] dx-t 2 J p t u ■ Vu tt ■ (L, Hfdx - t 2 j p t u • V(L, Hf 
= -t 2 j put-V [(L, Hf ■ u tt ] dx-t 2 j p t u ■ Vu tt ■ (L, Hfdx 

-t 2 J pu-V[u- V(L, Hf ■ u a ] dx := hi + hi + hs 

where the superscript ' means the transpose of the vector (L,H). 
Now, direct estimates yields that 



u tt dx 



\hi\ <t 2 \\^utth\\VPl 



,11V (L,Hf\\ 2 +t 2 



) iiv« tt || 2 ||ti t imi(L,fl) t 



< CV 



\\^u tt \\ 2 \\u t \\ H 2\\V{L,Hf\\ 2 + \\Vuu\\ 2 \\u t \\ m \\V{L,Hf\\ 2 



< at 2 \\Vu tt \\j + C a t 2 \\^p Utt \\ 2 2 \\V(L,Hf\\i + 



2 + 

H 2 + 



U\V(L,Hf\\ 2 



<ai'\\Vu t tf 2 + C a t^utt\\ 2 2 \\V(L,Hf\\ 2 + t z \\u t \\ 2 H2 + \\V(L,H) 



\t\\2 



where in the last inequality one has used Lemma |4,5 
Similarly, one can obtain 



I/22I <i 2 ||Vu tf || 2 |M 



<at 2 \\Vu tt \\l + C a r 



4(L,H)% 

2 m \\V{L,Hff 2 < at 2 \\Vu t t\\ 2 2 + C a \\V(L,Hf\\ 2 



and 



|/ 23 | < ^[lIvWIhllVHIocllVulUllVtL,^)*^ + \\Vuuh\\pu 2 \\oo\\V{L,Hf\\ 2 

+ llyP^I|2|IVP" 2 ||oo(||V 2 U i ||2 + |M|oo||V 3 u|| 2 + ||V«||oo||V 2 u|| 2 ) 

<at 2 \\Vutt\\ 2 2 + C a \\V(L,Hf\\ 2 



+c 



t'WVputt ||^(t||V«||^ + 1) + t z \\V 2 u t \\ z 2 + t z \\u\\ z H3 + ||V(L, Hf\\ 2 



Continuing, using the lemmas obtained so far, one can get that 



<t z 



Wy/puWooWVutthW^/pUtth + IIP^Iloo || Vn tt || 2 (||^i • Vtt|| 2 + \\u • Vu t || 2 ) 

+ \\VP u \\°o\\^u t th{\\V{ut ■ Vn)|| 2 + ||V(« • Vut)|| 2 ) 



llv^lloollVuttlWIv^ttlh + ||/Hloo||Vu tt ||; 



l||Vlt|| 4 + ||lt||oo||Vlti||; 



+ ||\/P"l| o||>/pUtt||2(KI|4||V 2 u||4 + ||Vlit||4||Vu||4 + |M|oo||V 2 U t | 



< ai z \\vu tt \\ z 2 + c a t z w^puuWi + ll^n^div^ni + 1) + 



2 

H 2 



(4.73) 



(4.74) 



(4.75) 



(4.76) 



35 



1-^4 1 < * 2 ||\/pwu||2||\/p||oo||«t||4||Vn t ||z 



(4.77) 



< Ct 2 \\^pu t t\\2\\u t \\ H 4V 2 u t \\ 2 < Ct 2 \\^pu tt \\ 2 2 + Ct 2 \\u t \\ 2 H1 \\V 2 u t \\l 

\h\ < C^HVputtlllllVulloo, (4-78) 



and 



|7 6 | <t ||divti tt || 2 ||A(p)t||4||Vn t ||4 < at \\divu u || 2 + C a t ||V u t \\ 2 , 
|J 7 | < t 2 ||divrt tt || 2 || A(p) t£ H2II Vw||oo < at 2 \\divu tt \\l + C a t 2 \\\(p) u \\l\\u\\ 2 H 3- 



(4.79) 
(4.80) 



Substituting the above estimates on Ii (i = 1,2,- •• ,7) into (|4.72j) and then integrating the 
resulting inequality with respect t over [r, t\] with r, t\ G [0, T] give that 



ti 



«?|lv^tt(<i)ll2 + y ^llVn^ll^^c + Cr 2 !!^^)^. (4.81) 
Since t^fputt G £ 2 ([0,T] x T 2 ) due to Lemma 1431 there exists a subsequence such that 



Tfc 0, 



rfcllvW^lll^O, 



as — >• +00. 



Letting r = T& in (|4.81|) and k — > +00, one gets that 



By, ([Q3]) . it holds that 



nlVpM*)ll2+ / s 2 \\Vu u (s)\\ 2 dt <C. 



sup [i 2 ||V 2 ^|| 2 (t)] < C sup t 2 ||v^«*t|l2(*)+* 2 lkllHi +* 2 IHIh3 + 1 

te[0,T] te[0,T] L 



< c. 



Finally, by (|4.66p . we can obtain 

sup [i 2 ||Vti|| 2 1/2 , 9 (i)] <C SU P 

tG[0,T] *G[0,T] 

So the proof of Lemma 14.61 is completed. 



t 2 \\u\\ 2 H 3 +t 2 \\u t \\ 2 H2 + 1 



< C. 



(4.82) 
(4.83) 
(4.84) 

(4.85) 
□ 



5 The proof of Theorem 11.1 



With the uniform-in-<5 bounds of the solution (p s ,u s ) in Lemmas 13. 1H3.7I and Lemma f4.lti4.61 one 
can prove the convergence of the sequence (p s ,u s ) to a limit (p,u) satisfying the same bounds 
as (p s ,u s ) as 5 tends to zero and the limit (p,u) is a unique solution to the original problem 
(|l.ip -( fL~4"j) , The details are omitted for brevity and one can refer to Cho-Kim [7] for the routine 
proofs. In the following, we will show that (p,u) satisfy the bounds in Theorem 11.11 and (p, u) 
is a classical solution to (Oj) . Since u G L 2 (0, T; H 3 (T 2 )) and u t G L 2 (0, T; 7T 1 (T 2 )), so the 
Sobolev's embedding theorem implies that 



u G C([0,T];H z (T z )) ^ C([0,T] x T ). 

Then it follows from (p,P(p)) G L°°(0, T; W 2 ' q (T 2 )) and (p,P(p)) t G L°°(0,T; ^(T 2 )) that 
(p, P(p)) G C([0, T];W 1,q (T 2 )) n C([0, T]; V4 /2,9 (T 2 ) - weak). This and then imply that 

(p,P(p))eC([0,T];W 2 ^(T 2 )). 
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Since for any r G (0, T), 

(Vu, V 2 u) G L°°(r, T; ^ X ' 9 (T 2 )), (Vu*, V 2 u t ) G L°°(r, T; L 2 (T 2 )). 

Therefore, 

(V«,V 2 m) G C([r,T] x T 2 ), 

Due to the fact that 

V(p,P(p)) G C([0,T];^(T 2 )) C([0,T] x T 2 ) 
and the continuity equation (|l.ip i . it holds that 

p t = u .\7p + pdivn G C([r, T] x T 2 ). 
It follows from the momentum equation (|l-ip 2 that 

(pu)t = C p u — div(pu ® u) — VP(p) 

= ij.Au + (p + A(p))V(divu) + (divtt)VA(p) + pu • Vu + pudivu + (it • Vp)u - VP(p) 
G C([t,T] x T 2 ). 
Thus we completed the proof of Theorem 11.11 



6 The proof of Theorem 11.2 



Based on Theorem ll.il one can prove Theorem 11.21 easily as follows. Since 

po G H 3 (T 2 ) ^ W 2 ' q {T 2 ) 

for any 2 < q < +00, it follows that under the conditions of Theorem II .2\ Theorem 11.11 holds for 
any 2 < q < +00. Thus, we need only to prove the higher order regularity presented in Theorem 
[L21 

Lemma 6.1 It holds that 
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sup 

te[o,T] L 



- / \\u\\ 2 Hi dt < C. 







VpV d u|| 2 (t) + ||(p,P(p),A(p))|| H 3(t) 

Proof: Applying d XjXk , j,k = 1,2, to (14.25D yields that 

pu XjXkt + pu ■ Vu XjXk + p X]Xk u t + p Xj u Xkt + p Xk u Xjt + p XjXk u ■ Vu + pu XjXk ■ Vu 

+Px 3 u Xk ■ Vu + p Xj u ■ Vu Xk + p Xk u Xj ■ Vu + p Xk u ■ Vu Xj + pu x . ■ Vu Xk + pu Xk ■ Vu Xj (6.1) 

+VP{p) X]Xk = pAu XjXk + V((p + X(p))divu) XjXk . 

Then multiplying (|6.1H by Au XjXk and integrating with respect to x over T 2 imply that 
J [p\^u XjXk \ 2 + V((p + \{p))&Yvu) X]Xk ■ Au X]Xk ]dx = j (pu XjXkt + pu ■ Vu XjXk ) ■ Au X]Xk dx 
+ J [P vk u t + p X] u Xkt + p Xk u Xjt + p XjXk u ■ Vu + pu X]Xk • Vu + p X] u Xk • Vu + p X] u • Vu Xk 
+Px k u X] ■ Vu + p Xk u ■ Vu X] + pu x . ■ Vu Xk + pu Xk ■ Vu Xj + VP(p) XjXk ] ■ Au XjXk dx. 

(6.2) 
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Integrations by part several times yield 

J (pu Xj x k t + pu ■ Vu XjXk ) ■ Au XjXk dx 



\Vu 



XjXfc I 



+ pu • V 



i=l 

^XiXjXfc ^ P^Xj ^XjXfc ^X^X 

i=l 

2 2 



i=l 



-Jt lp 



Nu XjXk \ 



-dx 



dx 



(6.3) 



[E 



+E. 

i=l 

PUxi ■ V //.,,,,. 



=1 i=l 
2 

da; 



iXjXfa 



i=i 



i=i 



and 



V((/i + X(p))divu) XjXk ■ Au X]Xk dx 

(p + A(/3))|V(divu) Xj:r J 2 + (X(p) XjXk V(divu) + X(p) x .V(divu) Xk + X(p) Xk V(divu) 



+ VX(p) x . Xk divu + VA(p)(d 



(6.4) 

Then substituting ()6.3j) and (j6.4j) into (|6.2p . summing over j, k = 1,2 and using the Cauchy and 
Young inequalities and the estimates in Sections 3-4, one has 



— I 

dt l 



^u\\ 2 2 + 2^||V^|| 2 (i) < C {\\uf H3 + l)\\(V i P(p),V d X(p))\\i + 1 



Next, applying d XtXjXk , i,j,k = 1,2, to (JH^ gives 

A^-^i + 3 a;iXjXfc (div(pn)) = 0. 



(6.5) 



(6.6) 



Multiplying ([63]) by p^^^,, 
to x over T 2 , one gets that 

d 



and summing over i,j,k = 1,2 and then integrating with respect 



dt 



W^pU < c\\v*p\\ 2 

<C||V 3 p|| 2 



||Vp||oo||V 3 «||2 + ||V 2 U|| 4 ||V 2 H| 4 + [|VT*||oo|[V 3 p|[2 + 

||V 3 n|| 2 + llVulUllVVlla + ||/o||oo||V 2 An||, 



>I|V 4 «||2 



<a||V^A^ + C a (||u||^3 + l)||V d HI 2 , 

where a > is a constant to be determined. 
Similarly, one can obtain 

d_ 

Let a = ^ . It follows from inequalities (|6.5p , (|6.7p and (|6.8p , that 



: ||(V 3 P(p),V 3 A(p))|| 2 <«||V 2 A«|| 2 + C Q (||n|| H 3 + l)||(V 3 P(p),V 3 A(p))||? 



(6.7) 



(6.8) 



|||(VpV 3 n,V 3 p,V 3 P(p),V 3 A(p))|| 2 (t) + p||V 2 An|| 2 (t) 

2 i 1 \ || /V73 _ vt3 t)( „\ w3 \ z' \ 1 1 2 



< c 



HI^ + l)||(V>,V d P(p),V d A(p))||^ + l 



(6.9) 
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Then integrating (|6.9p over [0, t) and using the Gronwall's inequality lead to that 



sup 

t£[0,T] 



\^pS7 6 uUt) + \\V 6 ( P ,P( P ),\(p))\\ 2 + \\V z Au\\ 2 dt<C. 



So the proof of Lemma [BTTI is completed. □ 

Now we prove other higher regularities presented in (jl.9p of Theorem ll.2i It follows easily 
from dSHD and that for any t x ,t 2 G [0,T], 



||VpV d «ll2(tl)-|l>/pV d «ll2(*2)-^0, 



(6.10) 



as ii — >■ i 2 - 

Thanks to Theoremll.il one has 



/9G C([0,T];H 2 (T 2 )) ^ C([0,T] x T 



(6.11) 



It holds that 



W\pV 6 u\mi) - \\pVM\t(h)\ = I / p'\V 6 u\ (h,x)dx- / //|VV(i 2 ,x)dx 



^^[^V 3 ^ 2 ^) - p\V 3 u\ 2 {t 2 ,x)]dx\ + | J p\V 3 u\ 2 {t 2 ,x) [p(h, x) - p(t 2 ,x)]dx\ 
< sup p(t,x) I / [p|V 3 u| 2 (ti,x) - p\V 3 u\ 2 (t 2 ,x)]dx\ 

+ sup / p\ V 3 u| 2 (t, x)dx sup \p(t\, x) — p(t 2 , x)\ 
te[o,T] J x&f 2 



p|V u\ (ti,x)dx— I p|V u| (t 2 ,x)dx + sup \p(t\,x) — p{t 2 ,x) 

xGT 2 



(6.12) 



< C 

— > 0, as t\ — > t 2 , 

where one has used ()6.10j) and ()6.11|) . 

Moreover, due to the facts that pV 3 u G L°°([0, T]; L 2 (T 2 )), p G C([0, T]; H 2 (T 2 )) and u G 
C([0,T];F 2 (T 2 )), it follows that pV 3 u G C([0,T];H 3 - w) which means that pV 3 u is weakly 
continuous with values in H 3 ((T 2 )). This, together with (I6.12D . leads to 



(6.13) 



(6.14) 



pV 3 u G C([0,T];L 2 (T 2 )). 

In a similar way, one can prove that 

(p,P(p)) G C([0,T];F 3 (T 2 )) ^ ^([O,^;^ 1 ^ 2 )). 

Moreover, since u G C([0, T]; # 2 (T 2 )) by Theorem O and p G C([0, T];H 3 (T 2 )) by (JEHD, one 
can prove that for any ii,i 2 G [0, T], 

||V 3 p U (ii,-)-V 3 p U (t 2 ,-)i^0, 

HVpVMt!,-)- VpV 2 n(t 2 ,-)||^0, 

\\X7 2 pX7u(t lr )-X7 2 pVu(t 2 ,-)\\ 2 2 ^0 



(6.15) 
(6.16) 
(6.17) 
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respectively as t\ — > t 2 . In fact, to prove (|6.17p . one has 
J \V 2 pVu{t u x) -V 2 pVu(t 2 ,x)\ 2 dx 

< J \V 2 p(ti,x) - V 2 p(t 2 ,x)\ 2 \Vu(t 1 ,x)\ 2 dx + J \V 2 p(t 2 ,x)\ 2 \Vu(ti,x) - Vu(t 2 ,x)\ 2 dx 

< C\\V 2 p(h, •) - V 2 p(t 2 , OHl + ||V 2 p||l||V«(ti, •) - V«(ta, -)ll4 

< C(||V 2 p(ti, •) - V 2 p(t 2 , -)||| + \\V 2 u(t u •) - V 2 u(t 2 , -> 0, 

as — > t 2 . Similarly, (|6.15|) and (|6.16|) can be proved. In view of (|6.13|) and (|6.15p - (|6.17|) . we 
have proved that 

pu G C([0,T]-H 3 {T 2 )). (6.18) 
The proof of Theorem 11.21 is completed. 
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